clear variables
syms C1 C2 R1 R2 U0 a
assume(0 < R1 & 0 < R2 & 0 < C1 & 0 < C2 & 0 < U0 & 0 < a)
Jednacine
syms C R
zamene = [R1 == R, R2 == R, C1 == C, C2 == C, a == 3]
zamene = 
vrednosti = [R == 1, C == 1, U0 == 10]
vrednosti = 
 
syms J U1 U2 V2 s
jednacine = subs([C1*(s*U1 - U0) + U1/R1 + (U1 - V2)/R2 == 0,...
-((U1 - V2)/R2) + C2*s*(V2 - U2) == 0, ...
J - C2*s*(V2 - U2) == 0, U2 == a*U1],...
lhs(zamene), rhs(zamene))
jednacine = 
 
promenljive = unique([U1, V2, U2, J])
promenljive = 
Resenje jednacina
odziv = solve(jednacine, promenljive)
odziv = struct with fields:
J: -(2*C^2*R*U0*s)/(C^2*R^2*s^2 + 1) U1: (C*R*U0*(C*R*s + 1))/(C^2*R^2*s^2 + 1) U2: (3*C*R*U0*(C*R*s + 1))/(C^2*R^2*s^2 + 1) V2: (C*R*U0*(3*C*R*s + 1))/(C^2*R^2*s^2 + 1)
 
odziv.U2
ans = 
syms u2(t)
u2(t) = ilaplace(odziv.U2, s, t)
u2(t) = 
simplify(u2(t))
ans = 
u2_num(t) = subs(u2(t), lhs(vrednosti), rhs(vrednosti))
u2_num(t) = 
Crtanje grafika
fplot(t, u2_num(t), [0, 10])
xlabel('t');
ylabel('u2(t)');
grid on;